Dissociative electron attachment to the H2O molecule. I. Complex-valued 
potential-energy surfaces for the ^Bi, '^Ai, and "^82 metastable states of the water anion 



Daniel J. Haxton,i'2>Q C. W. McCurdy,^'^ and T. N. Rescigno^ 

'Department of Chemistry, University of Califor-nia, Berkeley, California 94720 
'^Lawrence Berkeley National Laboratory, Chemical Sciences, Berkeley, California 94720 
''Departments of Applied Science and Chemistry, 
University of California, Davis, California 95616 

We present the results of calculations defining global, three-dimensional representations of the 
complex-valued potential-energy surfaces of the ^Bi, '^Ai, and metastable states of the water 
anion that underlie the physical process of dissociative electron attachment to water. The real 
part of the resonance energies is obtained from configuration-interaction calculations performed in a 
restricted Hilbert space, while the imaginary part of the energies (the widths) is derived from complex 
Kohn scattering calculations. A diabatization is performed on the ^Ai and ^B2 surfaces, due to the 
presence of a conical intersection between them. We discuss the implications that the shapes of the 
constructed potential-energy surfaces will have upon the nuclear dynamics of dissociative electron 
attachment to H2O. 

PACS numbers: 34.80.Ht 



I. INTRODUCTION 

Dissociative electron attachment (DEA) to the water 
molecule proceeds through a number of channels, each 
with a different energetic threshold, 



H,0 



H + OH- 
H2 + O- 

H- + OH {X 2n) 

H + H + 0- 
H- + OH* 
H- + H + 



3.27 eV 
3.56 eV 
4.35 eV 
8.04 eV 
8.38 eV 
8.75 eV 



(1) 



The production of these species occurs via three 
metastable Born-Oppenheimer electronic states of the 
H20~ system, whose vertical transition energies there- 
fore determine the incident energies at which DEA oc- 
curs. Those electronic states of the anion arc the ^Bi, 
'^Ai, and Feshbach resonances, and they are respon- 
sible for the three distinct peaks in the DEA cross section. 
Their potential-energy surfaces contain asymptotes cor- 
responding to the product channels listed in Eq.([T]), with 
the exception of the H-|-OH~ channel; this product is a 
result of nonadiabatic effects. 

Here we report the construction of the complex- valued 
adiabatic potential-energy surfaces associated with these 
resonance states, which will be used within the local com- 
plex potential (LCP) model [U, 0, S iSj] to calculate the 
nuclear dynamics leading to dissociation. The present 
study is followed by a second paper @, to which we will 
refer as paper H, in which we present the results of nu- 
clear dynamics calculations under the LCP model using 
the calculated surfaces. 
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Dissociative electron attachment to water was stud- 
ied as early as 1930, in the experiment of LozierQ, 
and as recently as 2006, in the study by Fcdor et 
aZ.Q. These two experiments, along with the rest of 
the prior theoretical and experimental work on this 

subject I, M M d, M, M [li, 13 113, M, H H IH, 

have succeeded in characterizing each of the product 
channels of Eq. ([T|) and the three Feshbach resonances in- 
volved in their production. However, prior to our recent 
theoretical study of DEA to water via the lowest-energy 
^Bi resonance [22! [23| . there had been no complete theo- 
retical treatment of this process, nor, in fact, any ab ini- 
tio treatment of dissociative attachment to any molecule, 
involving more than one nuclear degree of freedom. We 
will give a more complete summary of the prior theoreti- 
cal and experimental results concerning the dynamics of 
this process in paper H. 

The present treatment supersedes our previous study 
of DEA via the lowest-energy ^Bi state. We have also 
studied the angular dependence of DEA to H2O and 
H2S via the ^Bi state of either anion [23|. Addition- 
ally, we previously presented a qualitative study of the 
topology of the potential-energy surfaces of these three 
electronic states [2^, including the many intersections 
that these surfaces exhibit. That qualitative study in- 
forms the present study, in which we construct quantita- 
tive surfaces. 

Two of the most significant features of these surfaces 
can be seen in Fig. [1] reproduced here from Rcf. [25| . This 
figure depicts the behavior of the resonance energies with 
respect to the H-O-H bending angle 9hoh, fixing the OH 
bond lengths at ri=r2=1.81ao. The ^Bi and '^Ai reso- 
nances are degenerate at linear geometry; this degener- 
acy will lead to Renner- Teller coupling between the two 
states. In addition, there is a conical intersection be- 
tween the '^Ai and ^i?2 surface located at approximately 
&HOH=Ti° ■ This conical intersection plays a crucial role 
for the nuclear dynamics of DEA via the uppermost ^-62 
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FIG. 1: Real parts of resonance energies, in units of hartrees, 
for OH bond distance = l.Slao in C2v geometry, plotted with 
respect to bending angle, in degrees. 

resonance state. 

The presence of a conical intersection between the adi- 
abatic ^Ai and ^i?2 states means that the nuclear dy- 
namics calculations that employ these potential-energy 
surfaces must either explicitly include the coupling be- 
tween them, or be performed in a diabatic basis with a 
smooth coupling term. Wc perform an approximate dia- 
batization of the calculated adiabatic '^Ai and ^-62 (1 and 
2 "^A') states to produce diabatic ^Ai and ^£^2 surfaces. 
This diabatization is based upon the diagonalization of a 
particular symmetry operation described in Sec. IV C[ not 
the explicit minimization of derivative matrix elements. 

A resonance state may be characterized by a width, 
r, and an energy. En, which are functions of the inter- 
nuclear geometry q. These quantities define a complex 
potential surface V{q), 

V{q)^ER{q)-i^. (2) 

The width F is related to the lifetime as r = 1/F (we use 
atomic units throughout). For a triatomic, the internal 
coordinates q may be the set of bond-angle coordinates 
(''I, 0). 

We use separate techniques to define the two compo- 
nents of the potential-energy surfaces, Efj{q) and F((f). 
We employ the complex Kohn variational method (26l.[27l. 
m, [2^, [33, HH, m, m, m, 111] to perform scattering calcu- 
lations which include the effect of the electronic contin- 
uum upon the anion state. These calculations yield both 
En and F, but wc discard the value of Er obtained from 
them. In its place, we use bound-state configuration- 
interaction (CI) calculations to obtain Er. 

In the asymptotic regions, the resonances become 
bound electronic states, and therefore the use of bound- 
state methods is entirely appropriate in those regions of 
nuclear geometry. Near the Franck-Condon region where 
these states are resonances, our CI treatment restricts 



the included configuration space to eliminate the ground 
state electronic continuum from the calculation. Thus 
our CI calculations neglect the shift in Er due to cou- 
pling with that continuum, a well-known effect explained 
by the Fcshbach resonance formalism [s^ . Since the shift 
in Er by this coupling is generally of the same order as 
the width, F, this is an excellent approximation in re- 
gions where the resonance is narrow. The three resonance 
states we treat here have as their dominant configuration 
an electron attached to a singly excited configuration of 
the neutral target (Fcshbach resonances), and thus tend 
to have small widths. However, there are a some geome- 
tries where the width of one or the other of the two upper 
resonances is of the order of a few tenths of an eV, and 
in those regions the error in our calculated values of Er 
for that resonance are larger. 

We also define a ground-state neutral H2O potential 
surface. The resulting energetics of the anion surfaces 
relative to the ground state appear to reproduce the ver- 
tical transition energies and all but one of the two-body 
asymptotes of these resonances very well, without a rel- 
ative vertical adjustment between the neutral and an- 
ion calculations. However, as we will describe in detail, 
the CI calculations fail to produce one of the two-body 
breakup asymptotes on one resonance surface, and we 
are forced to employ an ad hoc patching procedure with 
another CI to correct the unphysical behavior. A similar 
patching procedure was necessary to correctly describe 
the three body breakup asymptotes, but dynamics lead- 
ing to them is not part of the study we will present in 
paper II. 

The data points calculated at a large set of nuclear 
geometries are assembled into global representations of 
the potential-energy surfaces. We construct the global 
representations of the real part, Er, and the width, F, 
separately. The global representations of the real parts of 
the adiabatic ^Bi and the diabatic ^Ai and surfaces 
are defined by a sum of an analytic fit and a spline of the 
residual error of this fit. We also define a global represen- 
tation of the off-diagonal coupling between the diabatic 
states. The global representations of Er for the adia- 
batic ^Ai and states are obtained as the eigenvalues 
of the 2x2 Hamiltonian matrix defined by the fitted dia- 
batic surfaces. The adiabatic-to-diabatic transformation 
angle is applied to the constructed width surfaces to ob- 
tain diabatic widths. The final result is a complete set of 
complex potential surfaces for the full dynamics calcula- 
tions we report in paper II. 

The outline of this paper is as follows. We begin in Sec. 
im with a description of the electronic structure of the 
resonance states and the basic features of their potential- 
energy surfaces, reviewin g th e results on the topology 
of those surfaces of Ref. [2^. In Sec. [TTIl we discuss 
the complex Kohn scattering calculations from which we 
obtain F(g), present the results of these calculations, and 
describe the construction of the global potential-energy 
surfaces from the individual data points. In Sec. |V] we 
do the same for the configuration-interaction calculations 
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defining Eniq). In Sec. IVII we present the constructed 
potential-energy surfaces and discuss their features. 



II. ELECTRONIC STRUCTURE OF THE 
RESONANCE STATES 

The electronic states that are primarily involved in dis- 
sociative electron attachment tqJii^P are the ^Bi,^Ai, 

ij, M M 

231 . At the 




and Bo Feshbach resonances 

H m, m 113, M, M, M, El, 

equilibrium geometry of neutral H2O, ri ~ r2 = l.Slao, 
= 104.5°, where oq is the Bohr radius 0.529189379 x 
10^^" m, these states are characterized by their dominant 
electronic configurations 

: la? 2al Ihl 3a? lh\ Aaf 

la? 2a? Ibl 3a} 1&? 4a? 

la? 2a? Ibl 3a? 16? 4a? 

such that each is described as [H20]~^4a?. The parent 
state of each corresponds to the analagous configuration 
[H20]~^4a}, which may be a singlet or a triplet. Since 
the resonance states are doublets, both the singlet and 
triplet states of the neutral may be considered the parent 
state to which the incident electron is attached. 

The peaks in the DEA cross section occur at approxi- 
mately 6.5, 8.4, and 11.8eV, respectively, at the equilib- 
rium geometry of the neutral; these values approximate 
the vertical transition energies of the resonances. The 
ionization energy of H2O is 12.621eV. The vertical tran- 
sition energy to excite the resonance is close this 
value and we cannot eliminate the possibility that it is 
physically above the ionization energy. It is possible, 
therefore, that the parent state of the ^^2 resonance is 
an autoionizing excited state of neutral II2O — at least at 
some geometries. 

From the equilibrium geometry of the neutral, the res- 
onance states may be followed adiabatically toward the 
two-body breakup arrangements in which the electron is 
attached to one of the fragments or H-OH or O-II2. We 
identified the proper asymptotes in these two arrange- 
ments previously [25|. In the first of these arrangements, 
the ^Bi and '^Ai states correlate with H~-|-OII (^11), 
while the ^B2 state correlates with H~-|-OH (^S), leav- 
ing the OH fragment excited. In the second, the ^Bi 
state correlates with ground-state H2-l-0~. The '^Ai 
state is found at a much higher energy, as H2 (triplet 
laglcru)+0~ , and does not have a bound two-body 
asymptote in this arrangement. 

An additional complication is that the ^B2 surface is 
inherently double- valued. As described in Ref. [2^, its 
lower asymptote is ground-state H2-l-0", while the upper 
asymptote of this surface in this arrangement is not elec- 
tronically bound, and corresponds to O ( \D) plus a reso- 
nant state of . In our previous study [25| , we did not 
determine which resonant state of is involved. The 
doublc-valuedncss of the '^B2 surface presents a problem 



for our treatment because we are unable to fully charac- 
terize both sheets of this surface for all geometries. We 
instead define a single surface which interpolates between 
the two sheets within the three-body breakup region. We 
will attempt to identify any discrepancies between the 
results here and the experimental ones which we might 
ascribe to this omission. However, it is likely that this 
feature of the adiabatic manifold only plays a significant 
role for three-body breakup, a channel which we will not 
consider in paper II. 

We do, however, include the ^Ai-^B2 conical intersec- 
tion in the current treatment, because it is critical for a 
description of the dynamics of dissociative electron at- 
tachment via the ^B2 state. This conical intersection is 
a consequence of the crossing of the I62 and 3ai orbital 
energies as the bond angle is varied, and has analogs in 
both the singlet and triplet Ai and B2 states of the neu- 
tral, which exhibit conical intersections near the present 
one. 

The symmetry labels ^Bi, ^Ai, and ^B2 are appro- 
priate when the H2O molecule has C2V symmetry, i.e., 
when the OH bond lengths are equal. When the OH 
bond lengths are unequal, the molecule belongs to the Cs 
point group. The appropriate symmetry labels in those 
geometries are '^A" , 1 '^A', and 2 ^A'. Due to the conical 
intersection, the lower (1) and upper (2) '^A' states may 
each correspond to either ^Ai or ^B2, depending on the 
bond angle. A complete discussion of the topology of the 
anion surfaces is given in Ref. psj . 



III. FIXED-NUCLEI ELECTRON SCATTERING 
CALCULATIONS 

The resonance positions and widths are extracted from 
the results of fixed-nuclei scattering calculations carried 
out at physical (real) energies. The scattering calcula- 
tions yield an S matrix whose energy dependence is ana- 
lyzed to determine the location of the resonance pole in 
the complex energy plane as a function of nuclear geom- 
etry. These calculations are fully ab initio. We use the 
cornplex Kohn variational method [2^, [13, H^, H^, HO, [3ll . 
[3^ . [ssl . [3^ . [3^ , which provides a stationary principle for 
the S matrix. Since detailed descriptions of the method 
have been given elsewhere [sj, HH , we will limit ourselves 
here to a brief summary to establish the terminology we 
will use to describe our numerical calculations. 

The (A^ + l)-electron scattering wavcfunction is repre- 
sented explicitly in this calculation, using the standard 
methods of quantum chemistry: one-electron molecu- 
lar orbitals are assembled as sums of multicenter con- 
tracted Cartesian Gaussian functions, products of which 
define iV-electron configuration state functions used to 
construct the target states. This basis, augmented with 
additional Gaussian functions and numerical continuum 
functions to describe the scattered waves, is also used to 
expand the {N + l)-electron wave function. 

Thus, key components of this calculation include the 
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appropriate choice of one-electron orbital and multiclcc- 
tron configuration bases defining the target states and 
the resonance state. Not only must the resonance state 
be accurately represented, but also the target states into 
which it decays. This requirement becomes more difficult 
with increasing resonance energy, as more target states 
become energetically accessible as decay channels. 



Representation of the wavefunction and matrix 
elements 



In our implementation of the complex Kohn variational 
method, the {N + l)-electron scattering wavefunction is 
expanded as 



(+) 



A 



^Xr(ri--riv)^r?](riv+i) 
r 



(3) 



The first sum in Eq. ^ is over target states xr(ri...rAf) 
explicitly included in a close-coupling expansion, which 
may be energetically open or closed. The antisym- 
mctrizer is denoted by A, and the scattered wave as- 
sociated with channel T is further expanded as 

lr> 



T, 



rr,, 

llomrrio 



G\l\kr,r)\/klr 



(4) 

for incoming boundary conditions in the target channel 
Tq with initial loruQ quantum numbers for the incident 
electron. The functions (pi{r) in Eq.(l4]), which we de- 
note as "scattering orbitals" , are Gaussian molecular or- 
bitals that are orthogonal to the "target orbitals" used 
to construct the target states xr- Thus, one of the first 
steps in the complex Kohn calculation is the partitioning 
of the one-electron Hilbert space into the sets of target 
and scattering orbitals. The functions Jim and in 
Eq. (jlj are constructed from products of radial functions 

(ji and f?!^-' ) times angular functions that are real- valued 
combinations of spherical harmonics Yim consistent with 
the spatial symmetry of the anion. These functions are 
then Schmidt orthogonalizcd to the target and scatter- 
ing molecular orbitals. The function ji is the regular 

Ricatti-Bessel function, while g^^'' is a numerically gen- 
erated continuum function that is regular at the origin 
and behaves asymptotically like the outgoing Riccatti- 
Hankel function: 



(5) 



It is obtained by solving the driven radial equation 

/ 02 



^ ^ + 5n '(^r,r) = rj;(fcrr) exp(-ar^), 

(6) 



subject to the stated boundary conditions. For all calcu- 
lations presented here, we use a = 0.04. 

For energetically closed channels, only the scattering 
orbitals (pi{r) are included in the sum in Eq.(|3]); the con- 
tinuum functions jim and gp|2j are not included. Thus, 
the calculation can give S matrices that are discontinu- 
ous across channel thresholds. To minimize this problem, 
we include additional diffuse Gaussian functions in the 
scattering orbital basis to represent the wavefunction in 
barely-bound channels. In any case, the resonance ener- 
gies can be formally discontinuous across channel thresh- 
olds; their discontinuities are not due to errors in the 
calculation, but instead are properties of the true reso- 
nance states (37} . 

The second sum in Eq.® is over square-integrable 
{N + l)-electron configurations 6^ constructed exclu- 
sively from target orbitals. For convenience we refer 
to the (TV -I- l)-electron configurations as the "Q space" 
and to the close-coupling part of the expansion of the 
wavefunction as the "F space" of the calculation. The 
P space may be further divided into the bound compo- 
nent — built from (N + l)-electron configurations incor- 
porating only Gaussian orbitals — and the "free" compo- 
nent, corresponding to the target states times continuum 
functions. 

Two approximations are made in calculating the ma- 
trix elements of the electronic Hamiltonian with respect 
to the antisymmetrized basis functions of Eq.Q. The 
exchange portion of all matrix elements within P space 
that involve "free" components, as well as the matrix el- 
ements between Q space and the free components of P 
space, are assumed to be zero. These approximations fol- 
low from the orthogonalization of the free functions to the 
bound molecular orbitals and the assumed completeness 
of the combined sets of target 4- scattering orbitals over 
the restricted region of space spanned by the target or- 
bitals. Errors associated with these approximations are 
minimized by keeping the target orbitals compact and 
by augmenting the set of scattering orbitals with func- 
tions that extend beyond the target orbitals so that the 
orthogonalizcd continuum functions arc separated from 
the target orbitals by a large region of space. Further 
details concerning these approximations can be found in 
Refs. (Hill. 



B. Target states and basis 

The description of a complex Kohn calculation requires 
the specification of the target states xr, the Q space 
configurations 0^, and, for the expansion of the channel 
eigenfunctions, the scattering orbitals ipi and the Im pairs 
included in the asymptotic partial wave expansion. We 
first turn our attention to the target orbitals and states. 

The one-electron orbital basis was constructed entirely 
from SCF orbitals. We began with the following primi- 
tive Gaussian basis. On the oxygen, we used Dunning's 
triple-C basis [s^ plus polarization (d) and Rydberg (s, p, 
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TABLE I: Fifteen lowest H2O target states for Kohn calculation: energies at the equilibrium geometry ri=r2=1.8ao, 9 = 104.5°, 
and coefficient of dominant configuration in CI expansion, compared to energies from van Harrevelt and van Hemert [sHi and 
from our previous complex Kohn study [2^. 



State 


Enerev (hartree) 

dii rrpTii" 




Excitation ener; 
Prpvf22l 


?V (eV) 

Ref f38l 

J- Lol . 


Dominant config. 
Config. Cocf. 
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0.9883 
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6.463 
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4ai 


0.9944 




-75 7Q0751 7 


6.829 


7.932 


7.63 


l&i 


4ai 


0.9960 




-75 79Q4SQ1 


8.496 


9.511 




l&i 


^ 262 


0.9966 


Mo 


7c; 79fifi944 
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Q fil 1 
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l&i 


^ 262 


0.9963 






8 7fi4 








— 4ai 


U.yooo 




7c; 71 1 498'^ 


8 qS7 






l&i 


^ 2fei 


0.9593 


Ml 


-75.7100565 


9.025 


10.534 


9.95 


I61 


^ 26i 


0.7224 


'B, 


-75.7072318 


9.101 






I61 


— > 5ai 


0.9959 


'B, 


-75.7060192 


9.134 






Ibi 


— > 5ai 


0.9946 


Ml 


-75.6876243 


9.635 






3ai 


— > 4ai 


0.7012 


^Bi 


-75.6863850 


9.669 






l&i 


— > 6ai 


0.9936 




-75.6826701 


9.770 






Ibi 


— > 6ai 


0.9951 




-75.6691870 


10.137 






l&i 


la2 


0.9952 




-75.6688136 


10.147 






l&i 


la2 


0.9940 



and d) functions, augmented with an s function with ex- 
ponent 0.0955 and ap function with exponent 0.774. On 
the hydrogen, we modified the double-C phis diffuse basis 
of Chipman [40|. Chipman's basis consists of four con- 
tracted s functions and one p function. The most diffuse 
s function has exponent 0.0483. We replaced the single 
p function with two p functions with exponents 0.55 and 
0.13, which were chosen to minimize the energy of the H~ 
anion given by full CI in this basis, which was -0.52190 
hartrees, corresponding to a hydrogen electron affinity of 
0.59591eV, which recovers most of the experimental value 
of 0.75419eV. The total size of the contracted target basis 
was 54. 

Using this basis, we first performed a two-shell, gen- 
eralized SCF calculation corresponding to the average of 
a ten-electron (neutral H2O) and nine-electron (cation 
H2O+) SCF calculation. The purpose of such a calcu- 
lation is to obtain a basis that can describe both the 
neutral and Feshbach resonance wavefunctions, recalling 
from the previous section that the Feshbach resonances 
may be approximately described as two 4ai electrons 
bound to different states of the cation core. Further de- 
tails about the generalized SCF calculation can be found 
in the EPAPS archive[41|. 

The set of target orbitals used to construct the neutral 
states xr included a total of eleven orbitals — the five SCF 
orbitals plus six of the virtual orbitals: three ai, labelled 
4ai - 6ai; the la2; the 2bi; and the 262. The target states 
were obtained from a restricted configuration-interaction 
calculation within this 11-orbital space, in which the lai 
and 2ai orbitals were constrained to be doubly occupied, 
and all single and double excitations from the I&2, 3ai, 
and I61 orbitals into the set of six virtual orbitals were 



included. This "all singles and doubles" calculation is 
designed to account for the single excitations describing 
the dominant configurations of the excited states, plus 
relaxation of the remaining orbitals to first order; it also 
can describe correlation in the target, particularly in the 
ground state. 

Near the equilibrium geometry of the neutral, this 
description of the target states puts the excited states 
slightly lower in energy, relative to the ground state, 
than appropriate, but otherwise represents key features 
of their potential-energy surfaces well. In Table |T1 we list 
the energies of the low-lying target states at the equilib- 
rium geometry of the neutral, along with the results of 
van Harrevelt and van Hemert (ssj . 

For the close-coupling expansion of the scattering 
wavefunctions, we selected sets of target states xr based 
on spin and spatial symmetry. At higher energies, a 
larger number of target states was required to converge 
the complex Kohn calculation, and therefore we used a 
different number of target states for calculations on the 
^A" i^Bi), the 1 ^A' i^Ai or 2^2), and the 2 ^A' {^82 
or "^Ai) states. For calculations of the '^A" resonance, we 
included 15 target states, which near equilibrium geome- 
try are the 15 lowest energy states listed in Table HI For 
the 1 ^A' state calculation, we added 14 additional states 
of A' symmetry, for a total of 29 states. For the 2 '^A' 
calculation, we added an additional eight A" states, for a 
total of 37 states. The calculations are all converged with 
respect to the target state expansion at the equilibrium 
geometry of the neutral, but there are some geometries 
at which the calculations are not fully converged. 

The (A^-l-l)-electron configurations 8^ are constructed 
from the target orbitals and describe correlation, relax- 
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ation, and the penetration of the incident electron into 
the target. For the present calculation, they include the 
dominant configurations which contribute to the resonant 
states. This set of configurations was again obtained by 
keeping the loi and 2ai orbitals doubly occupied and 
distributing the other seven electrons over the remaining 
target orbitals, including at least one but no more than 
three electrons in the (4ai, 5ai, 6ai, la2, 262) virtual 
space. 



C. Scattering orbitals and insertion basis 

Of the original 54 molecular orbitals, 11 were used for 
the description of the target and the remaining 43 or- 
bitals were included in the set of scattering orbitals ipi. 
To these 43 orbitals we added an additional set of con- 
tracted Gaussian functions, all centered on the oxygen. 
This " insertion basis" is included to improve the two ap- 
proximations mentioned above associated with the ne- 
glect of certain bound-free and free-free matrix elements. 

The insertion basis was constructed from even- 
tempered sets of eight p- and d-type primitive Gaussian 
functions, rji. (To avoid linear dependence, we do not 
include s-type Gaussians separately, but rather include 
all Cartesian components of the d-type functions.) The 
exponents started at 0.07 for the p-wave set and 0.075 
for the d-wave set, and in both cases the ratio of con- 
secutiveexponents was 0.8. From these sets of primitive 
functions, we constructed a set of six p and six d-type 
contracted, orthogonal functions that were designed to 
have the minimum overlap with the most diffuse func- 
tions in the target basis. To this end, we constructed 
a matrix representation, in the primitive insertion basis, 
of the operator P which projects onto the space spanned 
by the most diffuse functions included in the target basis, 
i.e.. 



Insertion basis 

Dunning Rydberg 
Bessel l<=0.5 




Pij = '^{Vi\(t>k){'l^k\r]j), 



(7) 



5 10 15 20 25 30 

radius from oxygen center (units of aQ) 

FIG. 2: Modufus squared of diffuse basis functions for 1=1 
(p waves) used in the Kohn calculations, centered on oxy- 
gen nucleus: sum of insertion basis (dashed fine), Rydberg 
p- function with exponent 0.028 (thiclc sofid fine), and orthog- 
onalized Bessef function jiikr) (thin solid fine). The Bessef 
function has fc=0.5aQ"^ and arbitrary normafization. See text 
for further expfanation. 



the continuum functions weil for energies up to about 12 
eV. In the current work, the exponents were chosen to 
provide as complete a basis as possible within approxi- 
mately 20 atomic units from the oxygen center. 

A pictorial representation of the basis set used in the 
Kohn calculations, which demonstrates the extent of the 
diffuse target, insertion, and continuum basis functions, 
is shown in Fig. [51 In this figure, the squared modulus of 
various members of the diffuse and continuum basis for 
l=\ (j) waves) is plotted with respect to the radial dis- 
tance from the oxygen nucleus. The most diffuse basis 
function used in the target calculations, with exponent 
0.028, is plotted along with the sum of the squared mod- 
uli of the insertion basis, ^'^d that of the 1=1 
Bessel function, orthogonalized to the insertion basis. 



where the (pk are the most diffuse functions of the tar- 
get basis, Schmidt-orthogonalizcd. These functions com- 
prised the two s- , two p- and one d-type target functions 
with the smallest exponents. 

Diagonalization of P gave two p- and two d-type vec- 
tors with eigenvalues close to unity, which were discarded. 
The other six p- and six d-type functions were re- 
tained as the insertion basis, giving a total of 6 x (3-1-6) 
= 54 additional Cartesian Gaussian functions and a total 
(target plus augmented scattering) set of 108 functions. 
The exponents and contraction coefficients of the inser- 
tion basis are listed in the EPAPS archive 41 1. 

We sfrould remark that the insertion basis we use is 
very similar to the exterior basis used by Nestmann and 
others [i^l for use in i?-matrix calculations. In the R- 
matrix method, the Gaussian basis must represent the 
outgoing wave within the finite i?-matrix box. Nest- 
mann's basis is designed for a 20-bohr box and represents 



D. Extracting resonance parameters 

The complex Kohn calculation was performed at 10 
to 20 energies around the resonance location, and the S 
matrices produced were fitted to a Breit-Wigner form. 



bg 



'^nlm'^n' V rn' 

E-E '' 



(8) 



R 



where the background S^^ is either linear or quadratic 
in E, to obtain the resonance energy Eji and the width 
F at each geometry, as well as the partial amplitudes 
7„i,„, labeled by the decay channel index nlm, where n 
is electronic target state and Im label the angular mo- 
mentum of the emitted electron. The modulus squared 
of the partial amplitudes corresponds to partial widths, 
and in particular, the partial width with respect to decay 
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into the electronic channel n is 

r„ = ^|7nim|^- (9) 

Im 

Unitarity of the 5" matrix implies r„ = F. 

Dissociative electron attachment to water occurs at in- 
cident electron energies sufficient to excite multiple states 
of the neutral H2O target. The definition of the reso- 
nance location and width in terms of the pole in the S 
matrix, as per Eq.®, poses a problem near target state 
thresholds. In our Born-Oppenheimer treatment of this 
process, we examine the behavior of the fixed-nuclei reso- 
nance width and location as the nuclear geometry is var- 
ied. At some nuclear geometries, the resonance location 
may cross an excited state threshold. 

According to analytic S matrix theory as described 
by Newton [sj, in such a situation it is generally not 
the same pole of the S matrix that is responsible for the 
resonance feature in the cross section both above and be- 
low a channel threshold. In such a case, in accordance 
with the formal theory, we observe that one pole quickly 
replaces another as the geometry is varied, so that the 
location of the pole identified with the resonance is effec- 
tively discontinuous near such a threshold. In contrast, 
the description of dissociative electron attachment under 
the Local Complex Potential modelfT], Q requires 

continuous potential-energy surfaces. As a result, we con- 
struct a global representation of the resonance width that 
smoothly interpolates through such discontinuities. 

IV. RESULTS OF SCATTERING 
CALCULATIONS 

At the equilibrium geometry of the neutral, the reso- 
nance positions of the ^Ai, and states were cal- 
culated to be 6.09, 8.41, and 11.97 eV, respectively, with 
widths of 10.31, 28.8, and 193 meV and partial widths 
with respect to decay to the ground state of 10.31, 10.30, 
and 9.135 meV. At equilibrium target geometry, each res- 
onance lies below its parent state. The ^Bi lies 370 meV 
below its ^Bi parent. The ^Ai lies 351 meV below its 
triplet parent, which in turn is only 223 meV below the 
next ^Ai state {Ih^^ 26}). (These two '^^1 states are 
in fact on the edge of an avoided crossing, at the equi- 
librium geometry of the neutral, and thus the binding 
energy of the Ai resonance versus its parent is slightly 
reduced by this avoided crossing.) The ^i?2 lies 1.831 eV 
above the first ^B2 state, which has a dominant config- 
uration 3a]^^262, and 471 meV below its triplet 162^^4aJ 
parent. The latter state is the 26th root of the target CI, 
and has a vertical energy of 12.438eV; this value is near 
the ionization potential of water, 12.621eV, which calls 
into question whether the ^i?2 parent is a true bound 
electronic state of the target at all, at the equilibrium 
geometry of the neutral. If not, it is likely that it exists 
instead as a low-lying Feshbach resonance in e~+Il20+ 
Coulomb scattering. 



We performed scattering calculations along 13 distinct 
one-dimensional cuts within nuclear configuration space. 
A dense set of points was chosen along each of these cuts 
to give an accurate and descriptive picture of the behav- 
ior of the resonances. The upper 2 ^A' resonance was 
examined along nine of the 13 cuts. Six of the 13 total 
cuts lie in C2V geometry: symmetric stretch (ri = r2) for 
eHOH=75°, 105°, and 150°; bend for ri=r2=1.81ao and 
2.41ao; and versus the Jacobi coordinate R for constant 
HH bond length r^f^f =1.40ao- The remaining cuts are 
bend at ri = 1.81, ri=2.41ao, and single-bond stretch at 
the six combinations of ri={ 1.81, 2.41ao } and Ohoh={ 
75, 105, 150° }. For the 2 ^A' surface we did not calculate 
the three cuts at 150°, nor the bending cut at (ri=:1.81, 
r2=2.41). 

We present results for the three cuts which pass 
through the equilibrium geometry of the neutral: as a 
fimction of bending angle 6 for ri=r2=1.81ao; as a func- 
tion of r2 for ri=:1.81ao and 6'=105°; and as a function of 
the equal bond lengths for symmetric stretch at 0=105°. 
Also we present data for the cut in Jacobi coordinates 
at r/fjj=1.4ao, C2V symmetry (7=7r/2), which docs not 
intersect the target equilibrium geometry. 

These calculations test the limits of our implementa- 
tion of the complex Kohn method, and therefore it is 
useful to have a measure of the performance of the cal- 
culation. One such measure is the modulus of the 5" ma- 
trices produced. The physical 5 matrix is unitary, with 
eigenphases having modulus 1 . At the equilibrium geom- 
etry of the neutral, at the center of the ^Bi resonance, 
the modulus of the most nonunitary calculated S matrix 
eigenphase was 1.005. 

The '^Ai calculation performed even better, yielding 
a value of 1.004 for the modulus of most nonunitary 
eigenphase. The corresponding value for the ^i?2 cal- 
culation, at higher energy with 24 open electronic states 
of the target and a 46x46 S matrix, was 1.2. When 
in C's symmetry the on-resonance ^i?2 calculation ob- 
tains a largest size of 120x120 channels at ri=1.61ao, 
r2=1.81ao, 0HOff = lO5°; at which point the S matrix is 
also significantly nonunitary. We note in passing that 
in all cases, the unitarity of the calculation is much im- 
proved for energies off-resonance. 



A. Bend, ri = r2 ~ l.Slao 

Plots of the resonance positions and target states en- 
ergies included in the Kohn calculation for this cut are 
shown in Fig. [31 The resonance energies are lines, and 
the target energies are connected dots. The large-dotted 
lines are those target states included in all calculations; 
the small-dotted lines are not included in the ^Bi calcu- 
lation (small dotted A' states) or not included in both 
the ^Bi and 1 ^A' calculation (small dotted A" states). 
Dots are filled for Ai and A2 symmetries, and open for 
Bi and B2 symmetries. 

Within this cut. the target state curves compare well 
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Resonances and A' singlets, r^=r2=1.81 (units of a^) 
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FIG. 3: Location of resonances (solid curves) and target states (dotted curves) at ri=r2=1.81 (units of ao), as function of 
bending angle Ohoh- Filled dots, Ai and A2 target states; empty dots, Bi and B2 target states. The dotted curves, marked 
by an asterisk, refer to a second ^Ai resonance discussed in the text. 



with the high-quality CI results of Harrevelt and van 
Hemertjl^ . As previously mentioned and shown in Table 
m at the equilibrium geometry of the neutral the Kohn 
target states are about leV below their proper location. 
However, aside from this shift, the shapes of these bend- 
ing potentials arc quite comperable to those of Ref . [11] . 
The 1 ^Bi curves from both calculations are quite simi- 
lar. Both calculations show an avoided crossing between 
the 2 "^Ai and 3 ^Ai states (near Ohoh = 105° in our cal- 
culations), where they change character between lbi^2bl 
and 3aj~"'^4a}, although the crossing is broader in our cal- 
culations. The 1 state undergoes several changes in 
character as it passes several avoided crossings, but it 
seems that the Kohn target curve is again simply shifted 
by about leV lower versus the results of Ref [33. The 
1 state is predominantly 3a^^2&2 at Ohoh = 180° 
and changes character to \h'^^2a\ at the sharp avoided 
crossing near 120°. The 1 state nearly parallels the 
1 ^Bi state, slightly less than 2cV above it, for both cal- 
culations. 

Along the cut shown in Fig. O the ^Bi Fcshbach res- 
onance stays below its '^Bi parent, except at very small 
bending angles Ohoh, and therefore is below all target 



states except the ground electronic state. The '^Ai state 
follows its triplet parent, '^Ai [B.20]3a^^4a\ configura- 
tion. At most angles this configuration describes well one 
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FIG. 4: Width of the ^Bi and 1 '^A' resonances (dots) with 
respect to bending angle 6 for ri=r2=1.81 (units of ag), with 
interpolated global representation (plain lines). 
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FIG. 5: Top, width of the 2 ^A' resonance (dots) with 
respect to bending angle 6 for ri=r2=1.81 (units of ao), 
with interpolated global representation (plain line). Bottom, 
zoom of width in ^Ai symmetry, with width of novel ^Ai 
([H20]lbJ'^26i4ai) Feshbach resonance. 



of the target states, except where there are avoided cross- 
ings. There are two avoided crossings involving this state: 
the one with the [H20]16j"^26}, around 9hoh^100° 
and 9.5eV, and also one with the [H20]16^^2&2 state 
near 9hoh=55° and 13.5eV. Near both of these cross- 
ings, there are discontinuities as the resonances cross the 
avoided state threshold. For the former avoided cross- 
ing, the discontinuity in the resonance position is not 
visible on the scale plotted in Fig. [31 but there is a sig- 
nificant discontinuity in the '^Ai width at Ohoh=Q5°7 
for the latter crossing, the discontinuity in position near 
(^HOH~50° is apparent as well. The resonance 
also follows its triplet parent configuration and parent 
state through avoided crossings with the [H20]laJ~^2fe2 
state (avoided crossing at =92°; discontinuity at 

dHOH=Q6°) and [H20]16j^^la2 state (avoided crossing at 
8hoh=7Q°', small discontinuity at 9hoh=83°). The for- 
mer discontinuity is quite large (~ 0.25eV). 

The widths for the ^Bi and 1 "^A' (^^2 or ^Ai) states 
along this cut are shown in Fig. [D Both the raw Kohn 
results (dotted lines) and an interpolated version (plain 



lines, defined later, in Sec. IIVE[) arc plotted. The 
width of the ^Bi state is relatively constant. The '^Ai 
state attains a relatively large width around 0=95° as 
it crosses multiple target states, but then the width de- 
creases before the '^Ai state intersects the ^i?2 state be- 
tween Ohoh^TO and 75°. At this geometry the widths 
of the '^Ai and ^B2 Feshbach resonances happen to be 
nearly equal, and so there is only a small discontinuity 
as the 1 ^j4' resonance changes symmetry from '^Ai to 
^B2 as the bond angle is decreased. In Fig. |4] adjacent 
data points are connected by line segments only if no tar- 
get states are crossed as the geometry is varied between 
them. Therefore one can see that at most crossings on 
the 1 ^A' surface the discontinuity in the resonance en- 
ergy is in fact small. 

The width of the 2 '^A' (^Ai or ^^2) resonance is plot- 
ted in Fig. [5] The ^i?2 state attains a large width 
(0.25eV, T=2.5fs) as it becomes the 2 ^A' resonance, and 
in particular as it rises above the B2 target states near 
Shoh=82° and 95°, at which points there are significant 
discontinuities in the location of the physical resonance 
pole. In particular, at 95° there is a large discontinu- 
ity in the real part of the resonance location: the real 
parts of the energies of the two poles avoid each other by 
approximately 0.3eV, whereas the widths avoid by only 
0.025eV. In Fig. [5] adjacent data points have been con- 
nected if the resonance does not cross a target state of 
the same spatial symmetry, or otherwise exhibit a large 
discontinuity. For the ^B2 Feshbach resonance, there are 
large discontinuities for non-i?2 crossings at the 2 Mi 
crossing near 6hoh^8^° , and at the crossing of both the 
5 ^Bi and 5 ^Bi near 61^0^=108°. 

The calculation supports a different Feshbach reso- 
nance of ^Ai symmetry for 9hoh < 80° along this 
cut, whose parent is the [H20]15j^^25i[ "^Ai target 
state, and which would have the dominant configura- 
tion [Il20]15j"^25}4a}. The location of this Feshbach 
resonance is also plotted in Fig. \3\ For these geome- 
tries it is impossible to say whether this is an artifact 
of the calculation or a physical state, since along this 
cut, this Feshbach resonance is extremely narrow and 
is not bound by more than 25meV. A small upward 
perturbation of the calculated resonance location would 
lead to its disappearance above its parent. This state is 
only present when the [H20]16j~^26i[ configuration is the 
lowest-energy '^Ai target configuration, for then it cannot 
decay to the triplet or singlet Bi ([Il20]16j^^4a}) target 
states. This is the case when the bond lengths are mod- 
est and the bond angle is small. The Jacobi coordinate 
cut, which we present in Sec. IIV D[ gives a better view 
of the interaction between this Feshbach resonance and 
the other '^Ai resonance and demonstrates that in fact 
this [H20]16j"^25}4a} Feshbach resonance may exist as 
an asymptote of the [Il20]3aj"^4af Feshbach resonance, 
due to branching of the adiabatic PES. 
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FIG. 6: Location of resonances (solid curves) and target states (dotted curves) at 6=105°, as function of symmetric stretch 
ri=r2- Filled dots, Ai and A2 target states; empty dots, Bi and B2. 



B. Symmetric stretch, ri = V2,0 = 105° 

The target energies and resonance positions along this 
cut are shown in Fig. [6l with the same conventions as 
Fig. [3l The (undotted) curves near the top right part 
of each panel correspond to the first target state of each 
symmetry that is not included in any of the Kohn calcula- 
tions. The l^Bi curve is very similar to that of Ref. [ssl ]. 
except that the minimum is at 1.95ao, not 2. lap. The 
1 ^ A2 asymptote of the Kohn target curve is too high by 
-1 eV. The 2 ^Ai and 1 ^B^ state curves are similar to 
those of Ref. [Ill, but shifted down ~0.75 eV, and in fact 
duplicate the apparent change of character of the 1 ^-62 
state at about ri—r2=2AaQ. 

As the symmetric stretch coordinate 7'i=7'2 is increased 
from the equilibrium geometry at LSlag, the ^Bi and 
'^Ai (1 "^A') resonances approach the bound 0+H+H~ 
asymptote of the three-body system and both become 
bound at approximately ri=r2=3.5ao. The '^B2 (2 '^A') 
resonance exhibits large discontinuities, due to both tar- 
get state crossings and the interaction with the ^i?2 shape 
resonance as discussed in Ref. [25|. At ri=r2>2.4ao, the 
lower-energy branch of the ^B2 shape-Feshbach system 



is uncovered beneath the 1 ^A2 state of the target, which 
correlates to 0+H+H. 

The widths for these cuts arc shown in Figs. [7] and 
[HI with the same conventions as in Figs. |4] and [S] For 



0.05 r 
0.04 




I • • -:t---^. , .„. 1 

1.5 2 2.5 3 3.5 

r.|=r2 (units of ag) 

FIG. 7: Width of the ^Bi and ^Ai (1 M') resonances (dots) 
with respect to symmetric stretch r\=r2, at 6=105°, with 
interpolated global representation (plain line). 
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FIG. 8: Width of the (2 ^A') resonance (dots) with re- 
spect to symmetric stretch, ri=r2, at ^=105°, with interpo- 
lated global representation (plain line). 

the 2^ A' state, the final branch of the state was not 
included in the global fit, and thus the interpolated value 
levels off near 0.2eV, corresponding to the width of the 
penultimate branch at ri=r2=2.6 ap. 



C. Single bond stretch, ri = l.Slao,^ = 105° 

The resonance locations and Kohn target state ener- 
gies along this cut are shown in Fig. [Q] Along this cut 
towards increasing r2, the 1 '^A" , 1 ^A", and 1 ^A' Kohn 
target states, which correspond to the 1 ^Bi, 1 ^Bi, and 

1 '^^1 states at the equilibrium geometry of the neutral, 
move downward in energy as they approach the ground- 
state H-f OH asymptote of the fragments. The 2 ^A' (2 
^^i) state surface is repelled by the ground state and 
is essentially constant, starting at 9 eV, with a gentle 
maximum at 9.25 cV at r2=2.5ao, and approaching 8 eV 
at the OH(^E)-t-H asymptote. The difference between 
this asymptote and the OH(^n)+H ground state of the 
fragments is 4eV for the Kohn states; this may be com- 
pared with an earlier calculated value [i^ of 5.27eV at 
ri = 1.80ao. The results for the excited state sing lets are 
consistently ^leV below the results of Ref. [Saj- Thus, 
the ground X ^Ai Kohn target state curve dissociation 
energy is underestimated by about leV along this cut. 

Moving towards increasing r2 along this cut, the ^A" 
and 1 ^A' Feshbach resonances fall in energy and be- 
come bound as H~-|-OH(X ^11) near r2=4.5ao. The 2 
^A' state correlates to the H~+OH(E) asymptote, which 
lies above the ground state of the neutral, H+OH(X ^H); 
its width therefore goes to zero only asymptotically. The 

2 ^A' state has two large discontinuities along this cut. 
The first, near r2=2.15ao, is associated with an avoided 
crossing between the '^A' parent state and another state 
that becomes the new parent, correlating to OH(^S)-|-H, 
triplet coupled. As r2 is increased further, the binding 
energy of the Feshbach resonance with respect to this 



state decreases to zero at r2=3.0ao, and it disappears 
for r2>3.0ao, but by this time a distinct resonance pole 
which follows the corresponding singlet state [correlating 
to OH(^E)-|-H, singlet coupled] has already appeared to 
take its place. 

Along this cut, the binding energies of the Feshbach 
resonances (which are plotted in the EPAPS archive) 
with respect to their parents all reach a minimum around 
r2=3.25ao. This is also the geometry at which the mid- 
dle 2 ^A' branch disappears. It is interesting that this is 
the case not only for the electronically similar ^A" (^Bi) 
and 1 '^A' (^Ai) states but also for the final branch of 
the 2 ^A' state. In each case this minimum is approxi- 
mately 50meV, which for the 1 and 2 ^A' states is only 
twice their width, F « 25meV. We certainly cannot be 
confident that our treatment of the TV- and (A^-l-l)- elec- 
tron systems is balanced at the 50meV level. This value 
is an order of magnitude smaller than the binding ener- 
gies at the equilibrium geometry of the neutral, and calls 
into question whether the physical states in fact may rise 
above their parents near r2=3.25ao. 

The widths for these cuts are shown in Figs. [10] and 
[TTl with the same conventions as in Figs. [3] and [51 The 
width of the ^A" state goes smoothly to zero as the res- 
onance becomes bound. In contrast, the 1 ^A' state 
achieves a large width as r2 increases, as high as O.lScV 
at r2=4.55ao. It may do so because it is connected to the 
ground state by an s-wave matrix element; we suspect 
that it does do so because at intermediate geometries, 
the electronic structure is highly correlated, and the or- 
bital description of the resonance is likely to be different 
from that of the target states, muddling the distinction 
between Feshbach and shape resonances. In other words, 
as the 1 ^A' Feshbach resonance dissociates it takes on 
an increasing degree of shape resonance character, lead- 
ing to an increase in width. We discuss this issue further 
in paper II. 



D. Jacobi coordinates, thh = 1.4ao,7 = 90° 

One may define a Jacobi coordinate system for II2O 
in which thh is the H— H bond length, R is the dis- 
tance between the oxygen and the H2 center of mass, 
and 7 is the angle between the two corresponding vectors. 
Here we present data along the cut r^f//=1.40ao, 7=90° 
{C2V geometry), plotted as a function of R. This cut 
traces the resonances from a squeezed geometry (small 
Ohoh) towards the 0-fH2 arrangement in C2V symme- 
try. It is uninteresting except for the interaction between 
the '^Ai (2 '^A') Feshbach resonance with configuration 
[H20]3aj~^4a^ and the second ^Ai Feshbach resonance 
mentioned above, with configuration [H20]16j'^26j4a}. 

In Fig. 1121 we show complex Kohn S matrix elements 
and fitted resonance locations along this cut. This figure 
requires some explanation. The target state energies are 
plotted as dots, and the resonance location is plotted as 
a single, thick solid line. At 15 of the plotted values of R 
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FIG. 9: Location of resonances (solid curves) and target states (dotted curves) at ri=1.81ao, 5=105°, as function of r2. 



we have also plotted the real part of the s ^ s S matrix 
element calculated using the complex Kohn method but 
with ordinate and abscissa reversed. The scale and ori- 
gin in the horizontal (ordinate) direction is arbitrary for 
these S matrix elements. 



The leftmost data in this figure are calculated at 
_R=1.67ao, a geometry that in bond-angle coordinates is 
approximately (ri=r2=1.81ao, 0ffo/f=45°). At this ge- 
ometry, both the [H20]3a^^4a2 Feshbach resonance, at 
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FIG. 10: Width of the ^A" CBi) (x 10) and 1 ^A' reso- 
nances (dots) with respect to r2 for ri = 1.81ao, 6'=105°, with 
interpolated global representation (plain line). 
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FIG. 11: Width of the 2 ^A' resonance (dots) with respect to 
r2 for ri=1.81ao, 5=105°, with interpolated global represen- 
tation (plain line). 
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FIG. 12: Location of Mi resonances and raw complex Kohn 
data at rHii=l-4ao, C2v symmetry (7=90°) as functions of 
Jacobi coordinate R. Solid thick line: resonance location. 
Filled dots: Ai target states. Empty dots: other target states. 
Thin lines: real part of s — > s 5 matrix element, on an arbi- 
trary scale, as a function of energy at various R, plotted with 
ordinate and abscissa reversed. 



~14.9eV, and the [YL20]lh^^2h\Aa\ Feshbach resonance, 
at ~12.75eV, are present and their Breit-Wigner profiles 
may be seen in the behavior of the S matrix element. As 
the coordinate R is increased beyond 1.87ao, however, 
the main [H2 0]3a]^^4a^ resonance is lost within a web of 
avoided and actual target state crossings around 14eV. 
The lower [B.20\lh~Hh\Aa\ Feshbach resonance is still 
present, however, and is bound by 0.5eV with respect to 
its [H20]16j"^26j ^Ax parent state at i?=1.97ao. This 
large binding energy brings up the possiblity that this 
resonance is in fact physical, at least at this geometry, 
and not simply an artifact of the calculation supported 
by recorrclation of the target. This Feshbach resonance 
follows the energy of its ^Ai parent through an avoided 
crossing near 2.1ao and through actual crossings with two 
^Ai states at approximately 2.45 and 2.55ao. 

At approximately i?=2.67ao there is an actual cross- 
ing between the '^^i parent and the ^Ai state, which 
is a parent of the familiar [H20]3aj~^4af Feshbach reso- 
nance, and at this point, the [H20]16^^26}4a} resonance 
appears to change character to that of the more familiar 
resonance, as there is a kink in the resonance trajectory 
and the resonance pole then follows the ^Ai parent as it 
dissociates toward H2 (lcrglcr„) -I- O ('^P). 

Thus, it appears that the "^Ai Feshbach resonance may 
undergo an interaction with a different Feshbach reso- 
nance of a sort similar to that which occurs psj within 
the B2 manifold between the [H20]l&2^^4aPB2 Feshbach 
resonance and a ^i32 shape resonance. The difference 
would be that in the current case, the topology is sup- 
ported by the disappearance and appearance of different 
branches of the adiabatic manifold due to the crossing of 
A^-electron target states, whereas in the ^i?2 case, it is 
supported by the underlying {N + l)-clectron Hamilto- 
nian. Another difference is that the topological complica- 
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FIG. 13: Width of the ^Bi and (1 ^A') resonances (dots) 
in C2v symmetry with respect to the Jacobi coordinate R 
for r//_f/=1.4ao, with interpolated global representation (plain 
line) . 



tions that may occur on the ^Ai potential-energy surface 
seem to do so at geometries not sampled by the propa- 
gating DEA wavcpacket. making them irrelevant to the 
physical problem, although we did not attempt a com- 
prehensive analysis of this issue. Along this cut, the ^Bi 
and ^-82 (1 "^A') states quickly become bound as 0~-|-H2, 
and we present the widths of these states in Fig. [T21 



E. Global representation of the widths 

For the purpose of performing nuclear dynamics calcu- 
lations, a global representation of the width F is required. 
Global representations were constructed separately for 
the ^A", 1 2yl', and 2 "^A' states. 

The first step in constructing each of these global rep- 
resentations was to define a continuous representation 
along each of the 13 one-dimensional cuts listed above, 
which was obtained via cubic splines in the coordinates 
ri, r2, and cos{6hoh)- Along these lines but beyond the 
last data points of the cuts in the ri, r2, or symmetric 
stretch directions, the width was either set to the termi- 
nal value, or to zero in the case of the large-r region for 
the asymptotically bound 1 ^A' or ^Bi resonances. 

The global representation of the fit was obtained as 
follows. Each one-dimensional cut is represented by a 
curve in the three dimensional space of coordinates q ~ 
(?'i,?'2,^)- For a desired geometry go, we first calculate 
for each cut (by interpolation, as explained above) the 
width at the point which is closest to go. We then take 
an average of the values at these 13 points, each inversely 
weighted according to its distance from q^. Examples of 
the interpolated witdh surfaces are shown in the EPAPS 
archive [41 1 . 



14 



V. CI CALCULATIONS FOR THE REAL PART 
OF RESONANCE SURFACES AND THE 
NEUTRAL SURFACE 

We construct the real part of the resonance energy Er 
as a function of the internal nuclear geometry of H2O 
using bound-state configuration-interaction calculations. 
Our task in doing so is to accurately reproduce real- 
valued component of the physical potential-energy sur- 
faces, taking into account the numerous features of these 
surfaces which were described in Ref. [25| . and, in par- 
ticular, the conical intersection between the and '^Ai 
surfaces, for which we require a diabatization. 

The main configuration-interaction calculations on the 
resonance states and the ground state of the neutral de- 
scribed below in Sees. IV Bl and IV Dl reproduce the ver- 
tical transition energies and all but one of the two-body 
breakup asymptotes correctly. However, the three-body 
asymptotes of the main configuration-interaction calcu- 
lation are all too high in energy (by as much as more 
than leV), and the asymptote of the diabatic '^Ai state 
is far below its proper 0^+H2(crgCr^) ^Y^u asymptote, 
instead being degenerate with the other resonances as 
0~ +ll2{(Jg)- For this reason, we must patch the ab ini- 
tio surfaces produced from the main CI calculation. 

We construct global representations of the adiabatic 
^Bi and the diabatic '^Ai and states, as well as 
the coupling, by fitting the data points produced from 
the main CI calculation at each nuclear geometry. The 
global representation of the adiabatic 1 ^A' surface is 
then defined as the lower eigenvalue obtained by diago- 
nalizing the 2x2 Hamiltonian matrix of the global rep- 
resentations of the diabatic surfaces and coupling. The 
1 ^A' surface will be used for calculations on DEA via 
the ^Ai (1 '^A') Feshbach resonance, because in that case 
the conical intersection is not expected to play a large 
role in the dyamics, and therefore the adiabatic basis is 
sufficient. (The corresponding vector potential is not in- 
cluded.) The diabatic surfaces and coupling are used to 
calculate cross sections for DEA via the higher-energy 
(2 ^A') state. 

The constructed global representations of the diabatic 
'^Ai and ^B2 and the adiabatic and 1 "^A' surfaces 
include the errors that we have already mentioned: the 
three-body asymptotes are too high, and the two-body 
asymptote of the diabatic ^Ai surface in the H24-0~ ar- 
rangement is too low. In order to fix these errors, we 
combined these global representations with other global 
surfaces designed to reproduce the correct behavior in the 
region in question. In each case a single patching surface 
is combined with a single CI surface such that the final 
surface reflects the correct behavior. The surfaces are 
combined by taking either the higher-energy of the two 
surfaces (for the ^Ai patching) or the lower-energy of the 
two surfaces (for the three-body patching), and smooth- 
ing the resulting cusps with a simple algebraic formula. 

The patching surface for the "^Ai diabatic state must 
correlate to the proper H2(lo'gl(T„)+0^ asymptote of 



this state. We perform another CI calculation that cor- 
relates to this asymptote, and construct a global patching 
surface from the results, using the combined analytic fit 
+ spline technique we employed for the resonance sur- 
faces. The patching surfaces for the three-body asymp- 
totes of the adiabatic ^Bi and 1 ^A' surfaces, as well as 
that of the diabatic ^-62 surface, take the form of a simple 
analytic potential. 

The diabatization of the 1 and 2 ^A' adiabatic CI roots 
is a requirement dictated by the nuclear dynamics calcu- 
lations. In the adiabatic basis, there are singular deriva- 
tive couplings near the conical intersection between the 
1 and 2 ^A' surfaces. We have not calculated these cou- 
plings from our CI wave functions, and in preparation for 
dynamics calculations on the coupled surfaces, we there- 
fore perform a diabatization upon the 1 and 2 ^A' CI 
roots to produce the set of diabatic ^Ai and ^B2 curves 
and the accompanying coupling potential. 

In Ref. [2^ , we described how a full characterization of 
the manifold of Feshbach resonances must also include a 
^i?2 shape resonance whose potential-energy curve inter- 
sects that of the ^B2 Feshbach resonance in branch-point 
fashion. We do not include this ^i?2 shape-Feshbach in- 
tersection, and instead define a single surface which in- 
terpolates between the two sheets within the three-body 
breakup region. Therefore, the nuclear dynamics on this 
surface is unlikely to accurately represent the dynam- 
ics leading to three-body breakup. The results that we 
will present in paper II indicate that three-body breakup 
probably comprises a large component of the cross sec- 
tion for dissociative attachment via the ^B2 state. 



A. Orbital basis 

We constructed a single orbital basis for all the CI cal- 
culations on the resonances. We began by augmenting 
the contracted Gaussian basis of Gil et al. [l^ with the 
following additional Gaussian functions: on the hydro- 
gens, two s functions with exponents 0.08 and 0.0333, 
and two p functions with exponents 0.2 and 0.05; on the 
oxygen, one s function with exponent 0.0316, and one p 
function with exponent 0.254. The basis comprised 77 
contracted Gaussians total. 

We first obtained an orbital basis by performing a 
symmetry-restricted SCF calculation on the ^Bi reso- 
nance, which yielded the lai, 2ai, I62, 3ai, l&i, and 4ai 
orbitals, labeled la', 2a', 3a', 4a', la", and 5a' in Cg 
symmetries. This SCF calculation is bound, i.e., it can- 
not decay to Il2 04-e~ by a symmetry-conserving rotation 
among the orbitals, because the l&i orbital is restricted 
to be singly occupied. The same statement is not true of 
the other resonances, which are described by a hole in an 
a' orbital, and therefore we were forced to use this SCF 
orbital basis for CI calculations on all three resonances, 
not just ^Bi. 

These SCF orbitals have avoided crossings (see the 
graph in the EPAPS archive [4l|). When the resonance is 
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fully dissociated, i.e., in the arrangement ~+H +H, the 
4ai (5a') SCF orbital is best described as ^Jlpl{'^ h^Is + 
^Hbis), i-C, the bonding combination of the two hydro- 
gen Is orbitals. These Is orbitals are more similar to the 
hydride Is orbital than to the hydrogen Is orbital. As the 
hydrogens are brought closer, the energy of this orbital 
decreases and eventually crosses the 3ai (4a') and 262 
(3a') orbitals. At finite 0~+H2 separation, these cross- 
ings are avoided. There are avoided crossings between 
the 3a' and 4a' (3ai and 4ai) orbitals near 9hoh=4:8° , 
and between the 4a' and 5a' (4ai and 262) orbitals near 
&HOH=(i^° ■ In C2V geometries (e.g., ri=r2=7.0ao) the 
crossing between the 4ai and 2&2 orbital is an actual 
crossing. 

In the definition of the multielectron configuration 
space that we use for the CI calculation, we treat the 
4ai (5a') orbital differently from the 3ai or 2^2 orbitals. 
Thus, it was useful to define a "diabatic" 4ai orbital 
whose energy smoothly crosses that of the others and is 
continuous in character. To this end, we performed a 
rotation among these three orbitals wherein we replaced 
the oxygen nucleus with an uncharged center, by obtain- 
ing the lowest eigenvalue of the Hamiltonian within 
the space of the 3ai, 4ai, and 262 orbitals. The ground 
state eigenvector of this calculation was defined as the 
diabatic 4ai orbital, and its complement as the diabatic 
3ai and 262 orbitals. The 4ai diabatic orbital is thus the 
"H2-like" orbital. Fortuitously, this rotation leaves the 
4ai orbital virtually unchanged at the equilibrium geom- 
etry of the neutral. In the OH-I-H" asymptote, the 4ai 
orbital is also left unchanged as the hydride Is. The ex- 
pectation of the Fock operator upon this diabatic 4ai or- 
bital (plotted in the EPAPS archive [iH) passes smoothly 
through the avoided crossings at 6hoh ~ 48° and 61°. 
We used this rotated SCF basis for CI calculations on 
the three resonances, and in the description of these cal- 
culations below, the notation "4ai" refers to the diabatic 
4ai orbital. 



B. Main CI calculation 

The main configuration-interaction (CI) calculation 
on the three Feshbach resonances employs a configura- 
tion space that is restricted to the configurations which 
contribute most to the resonances under study. This 
space was defined by the three resonance configurations 
[H20]16j;^4a^ , [H20]3aj;^4a? , and [HaOJlb^Ma?, plus 
all singles and doubles into the virtual space, with the 
lai orbital restricted to be doubly occupied. Thus, we ex- 
cluded from the configuration space all references wherein 
the combined occupancy of (I&2 3ai I61) was six. This 
restriction excludes configurations of the form [Il20]n^, 
which correspond to discretized continuum states of H2O 
-I- e~ . The total size of this configuration space is 111792 
in A' symmetry, and 106110 in A" symmetry. 

The restriction placed on the configuration space was 
critical in that it produced the resonance energies as low- 
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^A' roots of the main CI calculation at 



'"i=r2=1.81ao, as a function of bending angle 9. 



lying roots of the CI. The resonances were identified by 
their dominant configuration. The ^Bi (1 ^A") energy 
was found to be always the lowest root of its CI, and the 
'^Ai (1 ^A') state was almost always the lowest root of the 
CI of that symmetry. Therefore, the energies thereby ob- 
tained for these resonances were smoothly varying func- 
tions of the nuclear geometry, not suffering from avoided 
crossings with discretized continuum states. The energy 
of the ^i?2 state rises above discretized continuum states 
of the type [ll20]2bi^n^m^ and [H20]3aj~^n^m^ at cer- 
tain geometries, namely small OH bond lengths and near- 
linear H-O-H geometry. However, the resulting avoided 
crossings were generally observed to be very narrow. The 
discretized continuum states of excited [H20*]n^ are un- 
dercorrelated, relative to the resonance roots, by this 
multireference, all-doubles treatment, and therefore oc- 
cur at a higher energy than they would otherwise. 

In Fig. [131 we plot the ^A' roots of this configuration- 
interaction calculation as a function of bending angle 
Ohoh, at the equilibrium bond lengths ri=r2=1.81ao, 
showing the behavior of the CI roots near the conical 
intersection. Also apparent in this figure are the dis- 
cretized continuum states, many of which lie below the 
resonance root when the bond angle is large. 

The orbital and configuration basis for this CI calcu- 
lation is designed to describe well the resonances at the 
equilibrium geometry of the neutral and in the two-body 
dissociation channels. This calculation is not designed to 
reproduce the three-body asymptotes, and makes a large 
error in these regions. The 4ai and 262 orbitals have sig- 
nificant hydride Is character on both hydrogen centers, 
as opposed to hydrogen Is; since there is no relaxation of 
the hydride Is orbital included in the reference space, the 
double excitations into the virtual space must play that 
role, and the three-body asymptotes are therefore under- 
correlated and unphysically high in energy. We have not 
attempted a precise characterization of the CI roots in 
the three-body breakup region. 

We performed these CI calculations at various geome- 
tries on a grid based on H-H-0 bond angle coordinates 



16 



(not H-O-H) rHH: foH, and Oh ho- We used these coor- 
dinates so that the H2-l-0~ exit well (along with one of 
the OH+H^ exit wells) would be well-represented by the 
spline proceedure we use. In our previous treatment [2^ 
of the ^Bi resonance, we used H-O-H bond-angle coor- 
dinates, which led to an unphysically corrugated spline 
representation of the H2+0~ well. 

For this calculation, we defined a full grid of 23 r points 
between l.Ooo and 12.0ao and 29 9 points between 1° 
and 175°. This 23x23x29 grid includes 15341 points. 
The CI calculations were each performed at roughly 4200 
appropriately chosen points on this grid. 

The energies of the resonance CI roots at the equilib- 
rium geometry (ri=r2=1.81ao, 6'=105°) were -76.030888 
hartree (^Bi), -75.943508 hartree (Mi or 1 ^A'), and 
-75.802877 hartree (^Ss or 2 ^A'). 



C. Diabatization 

Because the ^Ai and ^i?2 states have a conical inter- 
section, and since we have not calculated the derivative 
couplings between them, a diabatization [4^ 1 is required 
for the nuclear dynamics calculations. Our method for 
performing this diabatization is an approximate method, 
based not on the explicit minimization of first-derivative 
matrix elements [i^ H^, IZtI IZsI . [49| but upon the diago- 
nalization of a property [50j to obtain smoothly behaved 
diabatic states. Our technique is thus analagous to a 
diabatization via the diagonalization of the dipole op- 
erator between states which undergo a charge-transfer 
avoided crossing [5l|, or the diagonahzation of the Iz 
angular momentum operator between adiabatic states 
which have a S-H conical intersection at linear nuclear 
geometry [H, [13] ■ 

The property we use for our diabatization is a symme- 
try operation — a reflection perpendicular to the molecu- 
lar plane — which is already diagonal in the adiabatic ba- 
sis both in C'2v geometries and in the asymptotic OH+H 
arrangement channel. In C2v geometries, this reflection 
is that defined by the plane that contains the C2v axis, 
and that is perpendicular to the molecular plane. In such 
geometries, the "^Ai state is an eigenfunction of this sym- 
metry operation with a eigenvalue of -1-1 with respect to 
this reflection, and the ^i?2 state has a eigenvalue -1. 
In the asymptotic OH-I-H arrangement, this reflection 
is defined by the plane perpendicular to the OH axis 
(which again is perpendicular to the molecular plane) 
that crosses through the oxygen nucleus. In these geome- 
tries, the matrix representation of the reflection operator 
in the basis of these states is diagonal, with the 1 ^A' 
state, which correlates to H~-|-OH(^H), having a posi- 
tive diagonal matrix element, and the 2 ^A' state, which 
correlates to H~-fOH(^S]), having a slightly negative di- 
agonal matrix element. 

In Fig. [13 we show the reflection plane for one C2v 
geometry and for one geometry approaching the H^+OH 
arrangement. The reflection plane is that plane that 




FIG. 15: Vectors involved in defining the reflection opera- 
tor whose diagonalization provides the '^Ai and ^B2 diabatic 
states from the 1 and 2 '^A' adiabatic states. The reflection 
plane is marked by a wide blue (grey) line and contains the 
vector V. 

contains the vector v and which is perpendicular to the 
molecular plane. The vector is a weighted sum of the 
unit vectors _Li and J-2- These vectors are defined as the 
unit vectors that are perpendicular to the corresponding 
OH bond vectors fl and r^, and that are contained in 
the HOH bond angle. The expression for v is 

V ^ exp(-ri/ro)-L2 + exp(-r2/ro)-Li, (10) 

where the length parameter ro is 1 bohr. This value 
was chosen on the basis of the dimensions of the current 
molecular system, and by inspecting the continuity of the 
diabatic potential-energy surfaces produced. 

D. CI calculation for the ground state potential 
surface 

The Born-Oppenhcimer potential-energy surface for 
the neutral molecule is required for the nuclear dynamics 
calculations, and to set the zero of energy for the disso- 
ciative attachment cross sections thereby produced. 

For the calculation of the neutral potential surface, 
we followed a prescription similar to that used in gen- 
erating the resonance surfaces. An SCF calculation on 
the neutral was performed, followed by a configuration- 
interaction calculation with all singles and doubles from 
the SCF configuration, keeping the lai orbital doubly 
occupied. The size of this CI calculation is 22215 config- 
urations in C2v symmetry. 

For the neutral, we used the physical H-O-H bond an- 
gle coordinates and a grid defined by roH={^Al, 1.61, 
1.81, 2.01, 2.21, 2.41, 2.61, 3.01, 3.41 gq}, and 9={60, 
75, 90, 105, 120, 135, 150, 165°}. The CI calculation was 
performed on each point on this grid. 
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The energy of the neutral CI calculation at the equilib- 
rium geometry (ri=r2=1.81ao, 0=105°) was -76.2900969 
hartree, yielding vertical excitation energies for the res- 
onances of 7.054, 9.431, and 13.258eV. (No zero-point 
energy is included.) 



E. Global representation 

To construct a global representation of these five CI 
potential-energy surfaces — the ^Bi, the diabatic ^Ai and 
^i?2, the patching ^Ai surface, and the ground state sur- 
face — a reference potential was first fit to the data, then 
subtracted from the computed points; the remainder was 
then fit with three-dimensional cubic splines. The sum 
of the reference fit plus splined remainder comprises the 
global fit, which coincides exactly with the calculated 
points. 

The functional forms of the reference potential Vres 
for the four resonance curves and Vneut for the neutral 
are given in the EPAPS archive[4l|. Root-mean-square 
errors of each fit were on the order of O.leV. 

The errors of each of the analytic fits were fit to cubic 
splines. Since we did not calculate the full grid of points, 
a multi-step splining proccedure was required. First, a 
series of one-dimensional splines, in the Ohho direction 
and then along the thh and roH directions, was per- 
formed to obtain the splined error at the remaining grid 
points. Second, the full grid of data thus constructed 
was fit to three-dimensional cubic splines, and added to 
the analytic fit to obtain the global representation. This 
proceedure yields the spline surface Vs- The global rep- 
resentation is Vres + Vs for the resonances or Vneut + Vs 
for the neutral. 



F. Representation of electronic coupling term and 
transformation of width to diabatic basis 

The electronic coupling matrix element between the 
^Ai and states was represented by a fit to a poly- 
nomial times Gaussian expansion in a rotated coordinate 
system. The explicit form can be found in the EPAPS 
archive [ilf. The RMS error of this fit was O.OSeV. 

This global fit of the coupling matrix element has a 
small remainder. Thus, while the diabatic '^Ai and 



surfaces pass exactly through the calculated points, the 
coupling surface, and thus the adiabatic surfaces ob- 
tained by a diagonalization of the electronic Hamiltonian 
thereby constructed, do not do so precisely. 



G. Patching of the surfaces 

The potential-energy surfaces constructed from the 
main CI calculation appear to reproduce the known fea- 
tures and all but one of the two-body asymptotes of the 
physical system, without recourse to an overall vertical 



adjustment in the relative position of the neutral an anion 
curves, or any other ad hoc adjustment. However, there 
are two regions in which the current configuration inter- 
action treatment fails to reproduce the physical energet- 
ics: for the high-energy il2{crgcru) '^Su 4- asymptote 
of the 2 ^A' (^Ai) surface, and in the three-body breakup 
region for all three resonances. 

Since the goal of the present study is to present the 
most physically accurate theoretical treatment of disso- 
ciative electron attachment within the local complex po- 
tential model, we correct these flaws in the surfaces by 
employing a patching procedure. In both cases a second 
surface with the desired characteristics is constructed and 
patched to the errant area. This patching is performed on 
the global fits Vres+Vs, not upon the original data points. 
The patching is performed by taking either the maxi- 
mum of the original and the patching surface (for the di- 
abatic ^Ai surface), or the minimum (for the three-body 
asymptotes), and smoothing the resultant cusps with a 
simple mathematical formula. This formula perserves the 
surfaces identically in the unpatchcd regions, and is de- 
scribed below. 



1. '^Ai patching surface 

All three roots of the main CI calculation correlate to 
H2-l-0~ in that arrangement channel. However, the cor- 
rect asymptote of the 2 ^A' state in that arrangement is 
H2(l(Tgl(T„)-f 0~ [2^. Therefore, we performed an addi- 
tional CI calculation that correlates to this state in that 
arrangment, and patched this surface to the diabatic '^Ai 
surface produced from the main CI, thereby correcting it. 
The diabatic ^^2 surface is left unchanged. 

The Gaussian basis and orbitals for this CI were ex- 
actly the same as those for the main CI, including the 
rotation of the 4ai orbital; the only difference was the 
choice of configurations. We included all single and dou- 
ble excitations from the configuration [H20]4ai[, keeping 
the Ifli orbital doubly occupied, and the 4ai orbital never 
doubly occupied. We took the lowest root of this CI. In 
the O+H2 arrangement, the (diabatic) 4ai orbital corre- 
lates to the H2 lag orbital, and therefore the lowest root 
is 0~-|-H2 (l(jgla„). Elsewhere, the lowest root corre- 
sponds to a discretized continuum state of H2 0-|-e~, and 
lies below the diabatic '^Ai surface. We performed this 
calculation on the same grid as the main CI calculation. 
Wc constructed a global representation of the patching 
surface employing our function Vres and a splined resid- 
ual Vs, just as we did for the resonance surfaces. 

This surface and the diabatic ^Ai surface intersect, and 
the upper surface contains the proper asymptotes for the 
physical '^Ai surface. It also contains cusps where the 
surfaces intersect. These cusps were smoothed by the 
following prescription. Given the separation AE and the 
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average E of these two surfaces, eV) defined as 

= \Eai - Epatchingl , , , , , , , 

(11) A = 0.5 + 12exp(-2rffH/3 + l), (12) 

-C^ — yi^Ai + ^patchi7ig)/ ^, 

and a geometry-dependent minimum separation A (in the patched '^Ai surface was obtained via 
I 



ma.x{EAi, Epatching) > J^^- 



2. Three-body asymptote patching 

In order to patch the three-body asymptotes of the 
constructed global representations, we defined the fol- 
lowing two functions of the OH bond lengths ri and r2'. 

Plow (x) = 8.4287 + 5.8739 exp [-1.4838 (x - 1.81)] 
Ph^gh (x) = 10.7140 + 3.5886 exp [-2.4287 {x - 1.81)] 
X = min(ri, r2), 

(14) 

in eV, which functions have the same gradient as the 
surface in the symmetric stretch coordinate but are 1.5eV 
above it at the equilibrium geometry of the neutral. In 
the three-body breakup region, these functions approach 
the values of their constant terms, which are chosen to 
represent the physical asymptotes of the system. 

The function piow is used to patch the asymptotes of 
the ^Bi and the 1 ^A' surfaces. The value of its constant 
term, 8.43eV, is in-between the asymptotes H~-|-IlH-0^ 
at 8.75eV and H-hH-hO" at 8.04eV. While, as discussed 
in Ref. [l^, the proper asymptote of these adiabatic 
electronic states is the lower of these, the higher may 
be reached by an excursion into the OH-I-II" two-body 
breakup channel that does not rigorously follow the adia- 
batic state in becoming H^+H-l-O. Therefore, we choose 
a compromise between these values to represent the phys- 
ical system. A more accurate treatment would use two 
surfaces, but we expect the present treatment to be suf- 
ficient for determining the two-body DEA cross sections. 

The function phigh is used to patch the diabatic 
^i?2 surface. Its asymptote is chosen as the energy of 
B'+R+Oi^D). Unlike the ^Bi and 1 '^A' surfaces, the 
(2 '^A') adiabatic surface is inherently double-valued 
within the inner regions of nuclear configuration space, 
and has two proper three-body asymptotes; H-|-H-|-0~ is 
the other asymptote. We have not attempted to charac- 
terize the full double-valued ^i?2 surface. The patching 
surface phigh corrects the ^i?2 three-body asymptote to 
the higher of its two physical values, in order to repro- 
duce the 0H(2E)-|-H~ two-body asymptote accurately. 
Therefore, the repulsive wall in the H2-l~0~ potential 



well extends above its physical value. It is hoped that 
this treatment does not alter the dynamics leading to the 
two-body dissociation channels to a significant degree. 

We combine the functions phigh and piow with the res- 
onance surfaces using the same equation, Eg. p^ . as 
we used for the '^Ai patching, with a constant value of 
A=1.0eV. 



H. Transformation of width to diabatic basis 

In order to transform the width surfaces produced from 
the scattering calculations, which are constructed in the 
adiabatic 1 ^A' and 2 ^A' basis, to the diabatic basis, 
the adiabatic-to-diabatic transformation matrix was con- 
structed using the patched diabatic ^Ai surface and the 
coupling calculated from the original diabatization of the 
unpatched surface. This adiabatic-to-diabatic transfor- 
mation matrix is therefore different from the transpose 
of the one that diabatized the adiabatic states from the 
main CI calculation. The diabatic width surfaces are de- 
fined as 



1 C 1 / \ 1 2 2A' / 



with the adiabatic-to-diabatic transformation matrix U 
expressed in terms of the angle 0' , 



U : 



cos(|-) sin(|-) 



(16) 



sin(^) cos(l^) J ' 
which is defined in terms of the patched surface V2^'^ as 

V2A' ~ V2Ji, 



cot{0') 



2C 



(17) 



This angle 6' is different than the original angle 9 that 
diagonalized the rcfiection operator. 



cot(6') 



2C 



(18) 
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FIG. 16: Comparison of scattering calculations and CI results along various cuts. Bold lines: results of scattering calculations; 
thin lines, CI results. 



where V2A1 is the original unpatched diabatic '^Ai sur- 
face. In particular, in the 0~-|-H2 asymptote, the patch- 
ing vastly increases the difference between the ^Ai and 
^i?2 surfaces, and thus the new adiabatic-to-diabatic 
transformation matrix is nearly unity there. 



VI. DESCRIPTION OF THE COMPLEX 
POTENTIAL SURFACES 

Several views of the diabatic surfaces can be found in 



the EPAPS archive [4 1|. Here we show one in the vicinity 
of the conical intersection, and then focus on the adia- 
batic surfaces only. 



I. Comparison with complex Kohn results 

The results of the CI calculations arc compared with 
the resonance locations obtained from the complex Kohn 
calculation, along various cuts, in Fig. 1161 In most cases 
this comparison is quite favorable, although some differ- 
ences are apparent. These results are presented in terms 
of the original, unpatched CI surfaces. 

The top three cuts, which each contain the equilib- 
rium geometry of the neutral, are in excellent agreement, 
except for the 2 ^A' surface for the symmetric stretch 
cut. Along this cut, the resonance has branched 
into the two components of the double- valued ^i?2 shape- 
Feshbach state, as discussed in Ref. psj . 

We include the cuts at ri=2.41ao in Fig. [TBI because 
these geometries are relevant to the wavepacket dynamics 
on the 2 '^A' surface. The gradient of the real component 
of that surface, as well as the behavior of the imaginary 
component, takes the propagated wavepacket through 
these cuts. The cut at ri=2.41ao, 0hoh=75° is near the 
conical intersection. Along this cut, the behavior of the 
CI surfaces mirrors the behavior of the Kohn resonances, 
though the agreement is not as good at 9hoh=W5° . 



A. Views of the conical intersection 

A plot of the conical intersection is shown in Fig. 
[T71 In this figure, the symmetric stretch coordinate is 
held constant at ri -I- r2=3.62ao; we see one cut of the 
potential-energy surface that intersects the conical inter- 
section seam at a point, at a bond angle of approximately 



B. Energetics of the adiabatic surfaces 

The real parts of the resonance surfaces Vres + Vs with 
patching are plotted in Fig. [THl fixing the bond lengths 
at their equilibrium values and varying the bond angle. 
The initial state starts at 9hoh~WA.5°, at which geom- 
etry the resonances are 6.63, 9.01, and 12.83eV above the 
ground vibrational state of the neutral. The conical inter- 
section is apparent at approximately 0hoh=76° , where 
the ^Ai and ^i?2 surfaces intersect. The 1 ^A' surface is 
plotted separately and is lower than the diabatic surfaces 
near the conical intersection. This is a consequence of the 
three-body patching procecdurc and the fact that differ- 
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FIG. 17: Cut of conical intersection at ri + r2=3.62ao: adi- 
abatic representation, top; diabatized surfaces and coupling, 
bottom. Top, values of 1 and 2 ^A' potential surfaces as 
constructed from global representations K-es + Vs of diabatic 
^Ai and surfaces and coupling. The potentials are plot- 
ted with respect to asymmetric stretch, units of bohr, and 
bending angle in degrees. The ^Ai and ^-62 surfaces along 
C2v geometry (ri=r2) are marked with bold lines. Bottom, 
fitted diabatized surfaces. The coupling is plotted as contours 
at bottom, contours every 0.25eV. 



ent functions Phigh and piow are used to patch the 1 ^A' 
and ^i?2 surfaces. The patching of the surface is ap- 
parent as the sUght downward kink of the surface near 
115°, at the edge of the Franck-Condon region of the 
neutral. The bump in the 1 ^A' and ^Ai surfaces near 
140° is a locahzed artifact of the sphne which thwarted 
removal. 

The vertical transition energy for the ^Bi state is very 
near the experimental peak maximum for DEA via this 
resonance, while the vertical transition energies for the 
'^Ai and ^i?2 states exceed the experimental DEA peak 
positions by '^0.4eV and '--^leV, respectively. We must 
point out, however, that the location of the experimen- 
tal peaks in the DEA cross section do not necessarily 
coincide with the vertical transition energies, especially 
for the upper states with their shorter lifetimes. As we 
will see in paper II, the vertical transition energies of 
the upper resonance states are probably closer to their 
appropriate physical values than than these comparisons 
would suggest. 

Globally, these potential-energy surfaces appear to re- 
produce the essential energetics of the underlying phys- 
ical states. Only in the case of the three-body breakup 
region is this agreement the result of an ad hoc procedure; 
elsewhere, the potential-energy surfaces represent the re- 
sults of ab initio calculations. In the case of the two- 
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FIG. 18: Potential-energy curves of calculated Feshbach res- 
onances at ri=r2=1.81ao, as constructed by global represen- 
tation. 



body asymptotes, we have been fortunate to obtain very 
good agreement with the proper energetics. The energet- 
ics of the system of three coupled Feshbach resonances 
as calculated is summarized in Fig. On the left of 
this figure are the accepted values for the differences in 
internal e nergy amon g e ach of these species, obtained 
from Refs.lialll, M,m, On the right are 

the results of the present calculations. The three-body 
asymptotes listed there correspond to the unpatched val- 
ues of the configuration-interaction surfaces at the geom- 
etry (ri,r2,6') = (10.0ao,10.0ao,60°). 



C. Asymptotes of the adiabatic surfaces 

The two-body asymptotes of the potential-energy sur- 
faces are plotted in Fig. [501 and compared with bench- 
mark theoretical calculations on the diatomic fragments. 
The dots in this figure represent the values of the global 
representations of the potential-energy surfaces evalu- 
ated along the cut in Jacobi coordinates at i?=10.0ao, 
7=90°. Also plotted in this figure are the values of bench- 
mark calculations for the diatomic H2 [6l| and OH (X 
^n) [g^I fragments, which are shifted so that their zero- 
point energies lie at the accepted energy above our calcu- 
lated ground vibrational state energy of neutral H2O; the 
solid horizontal lines represent the accepted energy of the 
three-body channels, again shifted to correspond to our 
ground state H2O energy. The theoretical calculations 
of Chu, Yoshiminc and Liu (60| slightly underestimate 
the true dissociation energy of OH (X ^H) and (^E). As 
is clear from the comparison with these benchmark cal- 
culations, the energetics of the two-body asymptotes of 
these anion surfaces are reproduced extremely well by 
our calculations. The ground vibrational state of each 
agrees with the accepted value to within O.OSeV. The 
three-body asymptotes have been adjusted by the patch- 
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FIG. 19: Left: physical thresholds [5^, 1^, [S^, l56|, |57|, |5g, ISg 

of one- and two-body breakup channels, relative to ground 
state neutral H2O, relevant to dissociative electron attach- 
ment to H2O. Vibrational ground states where applica- 
ble — zero point energies are included. Right: results of the 
present configuration interaction calculations on the anions, 
relative to the calculated ground vibrational state of the neu- 
tral. The vertical transition energies from the ground vibra- 
tional state of the neutral to each CI surface are marked with 
arrows, with experimental peak maxima in parenthesis for 
comparison. The three-body asymptotes labeled in the right 
panel correspond to values of global fits of potential-energy 
surfaces at (ri, r2, 0)=(lO.Oao, lO.Oao, 60°). 




ing procedure to correspond with the appropriate values. 



D. Complete views of the complex- valued adiabatic 
surfaces 

Complete views of the global fits of the real and imag- 
inary components of the adiabatic potential-energy sur- 
faces are shown in Figs. [^Tll^ The real and imaginary 
components are plotted together with contour lines, as 
a function of bond length, for various bond angles. For 
the imaginary component, the contour line nearest zero 
is bold, and subsequent contour lines depict the magni- 
tude of the imaginary component increasing. The con- 
tour lines for the real part, some of which are marked 
on the perimeter of each panel, correspond to the energy 
above the ground state H2O energy as calculated with 
our neutral H2O potential-energy surface, with the zero 
point energy included, and are therefore consistent with 

Fig. m 

In the view of the ^Bi surface in Fig. [^11 one can see 
the two OH+H^ channels at bottom right and upper left 
of each panel. In the bottom row of panels, which show 
the surface as the bond angle Ohoh is decreased, one can 
see the potential well which corresponds to the H2-I-O" 
channel develop along the symmetric stretch diagonal. 
This channel is the lowest energy asymptote; it reaches 



FIG. 20: Diatomic potential curve in OH - H asymptote (top) 
or H2 - O asymptote (bottom) of calculated Feshbach reso- 
nance curves at R=10.0ao (dots); benchmark theoretical cal- 
culations on OH (^H / ^E) [eOl and H2 [U (thin lines) and 
proper three-body asymptotes (thick lines), with zero point 
energy oriented relative to our calculated ground state H2O 
zero point energy by accepted values. 



below 4eV, as can be seen in the final panel at 15°. The 
bottoms of the H~-|-OH wells in the upper panels are 
above 4eV. Although the OH-|-H~ well does reach below 
4eV in the panels at 6'/foff =30° and 15°, it does so only 
at small OH-H separations, at which geometry there is a 
local minimum due to the dipole-anion interaction. 

The ^Bi resonance energy is relatively flat with bend- 
ing angle near the equilibrium geometry of the neutral, 
and this fact is apparent in the similar shape and value 
of the contour lines at small ri and r2 in Fig. [^T] from 
6'=150° to 61=75°. At 61=45° and beyond, the bending 
potential becomes repulsive and the contour lines move 
towards larger ri and r2. The flatness of the bending 
curve will cause the dissociating wavepacket to make min- 
imal excursions beyond the cut at 0=105°, and in par- 
ticular, only a very small fragment of the dissociating 
wavepacket will reach the 0~-|-H2 potential well at small 
9. The imaginary component of the ^Bi surface is small 
and localized near the Franck-Condon region. 
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FIG. 22: Complete view of 1 ^A' surface. Real part: solid contours, l.OeV spacing. Imaginary part (r/2): dashed contours, 
5meV spacing. The bold contour at 5meV is the lowest contour line for the imaginary part. Dot, 75°: intersection of conical 
intersection seam with this cut. 



The entire 1 ^A' surface is shown in Fig. [511 The 
real part of this surface is similar to that of the ^Bi 
surface, having the same asymptotes, but is everywhere 
higher (except at Hnear 9hoh=^&0° where they are de- 
generate), and is highest above the ^Bi surface along 
the conical intersection seam, that intersects this figure 
at 9hoh=75° at (ri=r2=2.3ao), which point is marked 



with a dot in the figure. The conical intersection seam is 
roughly parallel to the grain of the cut, and falls within 
the plane ri=r2. Thus, the dot marks the point where 
the 1 ^A' surface meets the 2 ^A' surface along this cut. 
Everywhere else, it is below. 

Near the Franck-Condon region of the neutral, the be- 
havior of the l^A' surface with respect to bending angle 
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FIG. 23: Complete view of 2 ^A' surface. Real part: solid contours, l.OeV spacing. Imaginary part (r/2): dashed contours, 
lOmeV spacing. The bold contour at lOmeV is the lowest contour line for the imaginary part. Dot, 75°: intersection of conical 
intersection seam with this cut. 



is different from the relatively flat ^Bi surface. In the cut 
at 0=150° these two surfaces are nearly degenerate, be- 
ing exactly degenerate at linear geometry. As the bend- 
ing angle is decreased, however, the 1 ^A' surface rises in 
energy, and the contour lines in Fig. [22] buckle outward; 
one such contour line reaches a point at the conical in- 
tersection at 9=75°. From this angle, the contour lines 
near the Franck-Condon region relax slightly downward 
going to 9=45°, which behavior traces the small well on 
the 1 ^A' surface apparent in Fig. [TBI around 9=55°. At 
smaller angles, the H2+0~ potential well develops. 

The imaginary component of this state is much differ- 
ent from that of the ^Bi state, being in general higher 
within the Franck-Condon region, and having the large 
(O.leV) peaks in the exit wells just as the resonance be- 
comes bound as H~-|-OH. These peaks are a consequence 
of the A' symmetry of the resonance, as discussed in 
Ref. [1^. They will lead to a large rate of autodetach- 
ment in this channel, and an isotope effect in the cross 
sections for H2O versus D2O that is disproportionate to 
the entrance amplitude, as discussed in Ref. [23|. They 
also may portend the breakdown of the LCP model near 
these geometries, for they indicate that the 1 ^A' state 
may exist as a virtual state there. We examine this issue 
more fully in paper II. 

The 2 ^A' surface is shown in Fig. [23] and is differ- 
ent from the other surfaces in every region. Its lowest 
points are along the well of the the H~-|-OH(^E) chan- 
nel where the surface drops to 8eV. This channel has 
a broader and shallower potential-energy well than the 
H~+OH(^n) asymptote of the other resonances, with 
a minimum at approximately roij=1.95ao. This sur- 



face exhibits a broad plateau in the (ri,r2) plane at 
&HOH=75° and lOeV. At this angle there is a dimple 
near the conical intersection at (ri=r2=2.3ao), which is 
again marked with a dot. There is no H2-l-0~ well evi- 
dent along the symmetric stretch direction at 9hoh^^5° 
because in this arrangement the 2 ^A' state is unbound 
as II2 (triplet l(Tgl(T„)-|-0^; II2 (triplet l(7glcr„) has a 
dissociative potential-energy curve. Thus, the panel at 
0HOH~^5° shows a cut along the top of this repulsive 
wall, going up past the last contour line at 16eV; this 
repulsive wall extends infinitely in the symmetric stretch 
direction at ever smaller 9hoh ■ Near the Franck-Condon 
region, one can see that the real part of the surface is 
high ('^15eV) at 9hoh=^5{)° , and slopes downward go- 
ing from panel to panel to 9hoh=75° , as the contour 
lines near ri=r2=1.81ao go from being convex to con- 
cave. Moving on to the panels at 45°, 30° and 15°, the 
repulsive wall of H2(lcrglcr,i)-f 0~ then quickly develops. 

This surface, constructed from the CI calculation, is 
single- valued, though the physical surface is not. The 
^i?2 shape resonance curve which intersects that of the 
^B2 Feshbach resonance in branch-point fashion[2^ has 
been omitted from the present treatment. The physi- 
cal 2 ^A' surface contains not only the Il2(lcrc,l(T„)-|-0~ 
asymptote (overall ^Ai symmetry) at small H-H bond 
distances, but also the metastable H^(1ct^ lau)+0{^D) 
asymptote (overall ^i?2 symmetry) as well; the two three- 
body asymptotes of this state are thus II+H^+0(^_D) 
and H-|-H-|-0~ (degenerate with 1 '^A'). We patch the 
three-body breakup region of the ^B2 surface to its up- 
per physical asymptote at lO.TleV. 

The imaginary component of the 2 ^A' state takes a 
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large value (~0.12eV) along the symmetric stretch direc- 
tion for 9hoh=^05° and 150° but otherwise is smaller; 
its representation is dominated by the interpolation be- 
tween the large-valued regions and the small-valued re- 
gions. In particular, in the H^-|-OH (^E) exit well the 
imaginary component drops to approximately 0.02eV by 
(ri = 1.8ao, r2=3.5ao). The interpolation is physical when 
ri or r2 is small. However, in the three-body channel, we 
interpolate between physically distinct sheets of the ^A' 
manifold. Thus, the imaginary component drops from 
lOOmeV near the diagonal for ri w r2 > S.Oao going from 
0hoh=W5° to 9hoh='!^° , but this behavior is unphys- 
ical, and corresponds to the interpolation between the 
two sheets of ^i?2 symmetry. The imaginary component 
remains small along the diagonal going from 0hoh=75° 
to 45° as the conical intersection is passed and the sym- 
metry of the 2^ A' state changes from to ^Ai. 



VII. CONCLUSION 

We have calculated the potential-energy surfaces nec- 
essary for a description of dissociative eletron attachment 



to H2O. In paper II, these surfaces are used in a study 
of the nuclear dynamics in the local complex potential 
model, and the cross sections for dissociative eletron at- 
tachment are calculated. 
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